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20.  Abstract^  (continued) 

/•''  r-y:  /  '  r  ^  •^■ 

This  first  portion  of  the  repor^escribes  two-component  velocity  measurements  and 
mean  concentration  measurements  in  the  region  where  a  concentrated,  drag-reducing, 
additive  solution  mixes  with  a  channel  flow  of  water.  The  additive  is  introduced  through 
flush  mounted,  angled  injectors  that  initially  place  the  solution  next  to  the  wall.  Near  the 
injectors,  the  mean  concentration  profiles  contain  very  sharp  gradients  that  separate  high 
concentration  wall  regions  from  a  low  concentration  core.  When  the  central  core  is  void 
of  additive  the  root^mean-square  of  the  normal  velocity  fluctuation,  V,  and  the  Reynolds 
shear  stress,  p  u^  are  unaffected  in  this  core  region.  These  two  quantities  are 
significantly  lower  in  the  high  concentration,  wail  region  compared  to  their  corresponding 
water  flow  levels.  Further  downstream,  where  the  additive  is  present  across  the  entire 
channel,  V  and  jp  uv  are  reduced  across  the  entire  channel.  This  indicates  that  once  the 
additive  has  imxed  beyond  the  wall  region  in  even  small  quantities  the  structure  of  the 
flow  is  altered  everywhere. 

The  second  portion  of  the  report  describes  single  component  laser  velocimeter 
measurements  in  a  Newtonian  channel  flow  where  the  Reynolds  number  was  50,000.  The 
purpose  of  these  measurements  was  to  deduce  the  average  time  between  bursts  at  this 
relatively  high  Reynolds  number,  ^fhe^  results  show  that  the  time  between  bursts  when 
scaled  with  wall  sliSar  velocity  and  kinematic  viscosity  has  a  constant  value  of  00  for  the 
Reyadds  number  range  9,400  to  50,000. 
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1.  INTRODUCTION 


The  basic  goals  of  this  project  are  to  determine  the  mechanisms  by  which  drag- 
reducing  additives  modify  the  turbulent  transport  near  walls  and  to  develop  optimum 
methods  for  injecting  these  additives  into  wall  bounded  flows  of  water.  The  purpose  is 
to  develop  methods  for  predicting,  controlling  and  manipulating  turbulent  wall  flows. 

The  addition  of  small  amounts  of  soluble,  high  molecular  weight  polymer 
molecules  to  water  flows  has  been  one  of  the  most  successful  methods  for  reducing 
drag  and  manipulating  the  near-wall,  turbulent  structure.  There  are  two  variations  in 
the  basic  technique.  In  the  original  method,  relatively  low  concentrations  of  the 
additive  were  injected  and  the  additive  mixed  at  the  molecular  level  with  the  water 
flow.  Recently,  Bewersdorfr  (1085),  Frings  (1085)  and  Berman  (1086)  have  reported 
results  for  flows  where  more  concentrated  solutions  were  injected  and  the  resulting 
flow  contained  one  large  or  many  small  threads  of  concentrated  polymer  solution. 
Currently,  this  project  is  concerned  only  with  the  original  technique  where  solutions 
with  concentrations  less  than  0.2%  are  injected  through  flush  mounted,  angled  wall 
slots  into  fully  developed  channel  flows  of  water  (see  Walker  et  al.,  1086). 

Two  somewhat  difiTerent  technical  issues  have  been  addressed  during  this  past 
year.  The  major  portion  of  this  report  will  describe  the  evolution  of  the  time-average 
concentration  proflles  and  velocity  statistics  downstream  of  the  optimal  injection 
process  described  by  Walker  ut  al.  (1086).  Additive  concentration  was  deduced  from 
an  optical  attenuation  technique  when  the  injected  material  was  marked  with  dye. 
Velocity  data  were  obtained  with  a  two-color,  two-component  laser  velocimeter  aligned 
to  measure  simultaneously  the  velocity  components  in  the  streamwise  and  wall  normal 
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directions.  These  data  describe  the  velocity  and  concentration  profiles  at  four,  eight, 
sixteen  and  thirty-two  channel  heights  downstream  of  the  injectors.  At  the  upstream 
position  the  additive  is  primarily  confined  to  a  region  near  the  wall.  At  the 
downstream  position  the  additive  is  well  mixed  throughout  the  channel  flow. 

The  other  technical  issue  centers  on  the  time  scale  associated  with  the  near-wall 
turbulent  structure  (the  time  between  bursts).  In  a  Newtonian  flow,  essentially  all  of 
the  turbulent  kinetic  energy  and  most  of  the  turbulent  transport  occurs  during  the 
burst  events  (see  Corino  and  Bradkey,  1069  and  Kim  et  al.,  1971).  The  burst  event  is 

I 

I  a  sudden  outrush  of  low  momentum  fluid  away  from  the  wall.  Associated  with  each 

I 

I  burst  is  a  sweep  or  inrush  of  high  momentum  fluid  toward  the  wall.  Since  the  burst 

event  is  associated  with  major  velocity  fluctuations  in  the  near-wall  region  where 
drag-reducing  solutions  have  an  effect  on  the  time-average  flow  field,  the  changes  that 
take  place  in  the  burst  events  in  drag-reducing  flows  are  of  particular  interest.  The 
burst  data  reported  here  were  acquired  at  a  Reynolds  Number  of  50,000  to  confirm  the 
inner  variable  scaling  discussed  by  Luchik  and  Tiederman  (1987)  for  Newtonian  flows 
and  are  needed  before  further  comparisons  with  drag-reducing  measurements  are 
made. 
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2.  MIXING  EXPERIMENTS 


These  initial  experiments  were  conducted  to  demonstrate  the  feasibility  of  the 
techniques  and  to  determine  the  evolution  of  the  mean  polymer  concentration  profile 
and  turbulent  velocity  statistics  downstream  of  a  wall  injection  of  a  polymer  solution. 
These  experiments  will  serve  as  a  guide  in  the  design  of  more  detailed  experiments 
that  will  follow. 


2.1  Apparatus  and  Procedures 


The  test  section  of  the  recirculating  water  flow  loop  used  in  these  experiments 
was  the  2.5  cm  high  by  25  cm  wide  rectangular  cross-section  channel  described  by 
Walker  et  al.  (1986).  The  channel  was  constructed  from  one  half  inch  acrylic  sheet 
and  is  more  than  one  hundred  channel  heights  long.  Polymer  solutions  are  injected 
through  flush-mounted,  angled  slots  located  in  both  long  walls  of  the  channel.  The 
average  slot  angle  is  25  *  with  respect  to  the  wall  and  the  streamwise  width  at  the 
wall  is  0.100  inches.  The  injectors  are  located  more  than  one  hundred  channel  heights 
downstream  of  the  inlet. 

The  flow  loop  is  driven  by  four  ninety  gallon  per  minute  centrifugal  pumps 
operating  in  parallel.  At  each  end  of  the  test  section  there  is  a  large  stilling  tank  to 
isolate  the  test  section  from  any  hydrodynamic  disturbances  in  the  flow  loop.  The 
upstream  stilling  tank  contains  a  perforated  plate  followed  by  screen  and  open-cell 
sponge  section  and  a  smooth  two-dimensional  contraction  at  the  outlet.  The  inlet  of 
the  channel  is  preceded  by  a  flow  straightener  consisting  of  closely  packed  plastic 
drinking  straws  which  insures  that  no  large-scale  vorticity  exists  in  the  channel  entry 


flow.  The  downstream  tank  contains  a  perforated  plate  to  damp  out  disturbances  and 
a  copper  coil  through  which  cooling  water  is  passed  to  maintain  the  channel  water  at 
a  constant  temperature. 

The  reduction  in  wall  shear  stress  which  occurs  during  polymer  injection,  was 
deduced  in  two  ways.  If  the  flow  is  approximately  fully  developed,  then  the  wall  shear 
stress  is  proportional  to  the  pressure  drop  and  the  fractional  reduction  in  pressure 
drop  resulting  from  polymer  injection  is  equal  to  the  drag  reduction.  Static  pressure 
drop  measurements  were  made  with  Gilmont  micromanometers  using  carbon 
tetrachloride  as  the  working  fluid.  Because  the  fluid  is  recirculated  with  intermittent 
injection  of  polymer  ,  the  water-flow  pressure  gradient  had  to  be  monitored  periodically 
during  an  experiment  to  insure  that  drag  reduction  due  to  polymer  accumulating  in 
the  flow  loop  did  not  occur.  These  checks  showed  that  once  the  polymer  solution  had 
passed  the  test  section,  it  was  no  longer  an  efiective  drag  reducer. 

The  second  method  of  deducing  drag  reduction  was  by  estimating  the  mean 
velocity  gradient  at  the  wall  from  near-wall  velocity  measurements  and  measuring  the 
mean  near-wall  polymer  concentration.  These  results  were  then  used,  along  with 
viscometric  data  to  determine  the  wall  shear  stress  directly  from 

Two  component  velocity  measurements  were  made  using  a  Thermo-Systems 
Incorporated  (TSI)  model  9100-6  laser  velocimeter  system  incorporating  a  500  mW 
Lexel  model  85  Argon-ion  laser.  This  two  channel,  two  color,  four  beam  system  was 
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oriented  so  that  the  optical  axis  was  parallel  to  the  long  walls  of  the  channel  and 
perpendicular  to  the  flow  direction  allowing  direct  measurement  of  the  streamwise  and 
normal  velocity  components.  To  allow  measurements  very  near  the  wall  the  lower 
beam  of  the  pair  used  to  measure  the  normal  velocity  was  displaced  to  the  optical 
axis.  The  beams  were  focused  using  a  250  mm  focal  length  lens  yielding  a 
measurement  volume  with  a  diameter  of  65  (xm  and  a  length  of  300  fxm. 

The  scattered  radiation  was  collected  in  forward  scatter  slightly  oET  the  optical 
axb.  Combining  a  250  nun  focal  length  lens  and  a  130  fxm  aperture  in  a  TSI  model 
9143  field  stop  assembly  results  in  a  field  of  view  about  300  /xm  in  diameter  for  the 
receiving  optics.  This  results  in  the  effective  probe  volume  length  of  about  1.1  mm  in 
the  spanwise  direction. 

Both  channels  of  the  laser  velocimeter  were  frequency-shifted  40  MHz  and  the 
photomultiplier  output  for  each  channel  was  electronically  down-mixed  to  yield  a  net 
frequency  shift  of  0.5  MHz.  Velocity  information  was  extracted  from  the  down-mixed 
photomultiplier  output  using  two  TSI  model  1980  counter  type  processors  operating  in 
the  N-cycle  mode.  Temporal  coincidenced  was  imposed  in  the  measurements  and  the 
data  were  transferred  to  a  Digital  Equipment  Corp.  PDF  11 /03  mini-computer  using  a 
TSI  model  1998  interface.  Data  were  then  transferred  to  a  VAX  11/780  computer  for 
statistical  analysis. 

In  these  experiments,  velocity  measurements  were  made  both  in  pure  water  flows 
and  in  flows  where  turbulent  fluctuations  in  the  polymer  concentration  occurred.  In 
order  to  eliminate  bias  in  the  velocity  statistics  caused  by  concentration  fluctuations, 
the  seed  particle  concentration  in  the  channel  water  had  to  match  that  in  the  injected 


polymer  solution.  To  accomplish  this,  the  channel  water  was  seeded  using  fat  particles 
from  whole  milk  so  that  the  particle  arrival  rate  did  not  change  significantly  when  the 
polymer  injectors  were  turned  on.  Ensembles  of  four  thousand  points,  acquired  by 
making  only  one  measurement  per  Doppler  burst,  were  used  to  calculate  velocity 
statistics.  The  two-dimensional  weighting  method  of  Tiederman  (1979)  was  used  to 
correct  these  ensembles  for  velocity  bias  since  the  sampling  rate  and  the  particle 
arrival  rate  were  correlated. 

Mean  polymer  concentration  profiles  were  measured  by  dyeing  the  injected 
polymer  solution  with  Rhodamine  B  dye  and  then  measuring  the  attenuation  of  a  laser 
beam  directed  across  the  span  of  the  channel  parallel  to  the  wall  when  the  injectors 
were  turned  on.  The  green  (514.5  nm)  line  of  a  Lexel  model  85  laser  was  used  and  the 
change  in  the  beam  intensity  was  measured  using  an  Oriel  model  71841  silicon 
photodiode.  This  combination  yielded  spatial  resolution  of  0.5  mm  in  the  direction 
normal  to  the  wall.  The  logarithm  of  the  intensity  change  is  proportional  to  the 
space-averaged  dye  concentration  along  the  path  of  the  beam.  Since  the  span  of  the 
channel  was  not  long  enough  to  result  in  a  temporally  constant  space-averaged 
polymer  concentration,  the  output  of  the  photodiode  was  input  to  a  logarithmic 
amplifier  to  obtain  a  signal  that  was  proportional  to  the  space-averaged  concentration 
along  the  beam  path.  This  signal  was  then  time  averaged  using  an  integrating 
voltmeter. 


Experiments  were  conducted  at  a  Reynolds  number  of  18500  based  on  mass- 
averaged  velocity  and  channel  height  in  the  2.5  cm  channel.  This  resulted  in  a  shear 
velocity,  u^  of  3.79  cm/s.  A  seven  hundred  ppm  aqueous  solution  of  SEPARAN  AP- 
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273  was  injected  through  the  flush-mounted  angled  injectors  at  a  flowrate  equal  to 
that  in  the  undisturbed  linear  sublayer  of  the  water  flow  (450  ml/min).  Mean  polymer 
concentration  measurements  were  made  at  locations  4,  8,  16  and  32  channel  heights 
downstream  of  the  injectors.  Two-component  velocity  measurements  were  made  at 
these  same  locations  both  with  and  without  polymer  injection. 

2.2  Results 

Figures  1  through  3  show  the  first  and  second  statistical  moments  of  the 
streamwise  and  normal  velocity  components  normalized  with  inner  variables  (shear 
velocity,  u^  =  'Sjr^jp  and  kinematic  viscosity,  u  for  the  flow  without  polymer 
injection)  at  the  four  measurement  locations.  These  measurements  indicate  that  the 
flow  was  typical  of  fully-developed,  two-dimensional  channel  flow. 

Variation  of  drag  reduction  with  position  downstream  of  the  injector  is  presented 
in  Figure  4.  These  data  were  calculated  from  both  pressure  drop  and  velocity 
measurements.  The  drag  reduction  levels  calculated  using  pressure  drop 
measurements  were  in  reasonable  agreement  with  those  obtained  by  Walker  et  al. 
(1986)  in  the  same  channel  at  a  Reynolds  number  of  17,800  and  the  same  injection 
flowrate  and  polymer  concentration.  Comparison  of  the  values  of  drag  reduction 
calculated  from  the  slope  of  measured  mean  velocity  profiles  at  the  wall  and  pressure 
drop  measurements  indicates  that  the  pressure  drop  measurements  give  a  reasonable 
approximation  to  the  drag-reduction  levels  calculated  from  the  average  wall  shear 


rate. 
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Figure  1. 


Mean  velocity  profiles  for  water  flow;  □,  x/h  =  4;  0»  8  ;  A, 
16;  V,  32. 


2. 


Root  -  mean  -  square  streamwise  and  normal  velocity 
profiles  for  water  flow;  □,  x/h  =  4;  0>  ^2. 
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Figure  3.  Turbulent  shear  stress  profiles  for  water  fiow;  □,  : 

O,  8;  A,  16;  V,  32; - ,  1  -  2  y/h  -  dU'^/dy'^. 
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Mean  polymer  concentration  profiles  are  shown  in  Figure  5  for  the  four 
measurement  locations.  At  the  near-injector  location  (x/h  =  4),  the  polymer 
concentration  near  the  wall  has  been  reduced  to  less  than  one-third  the  injected 
polymer  concentration  but  there  is  still  a  sharp  gradient  between  the  high 
concentration  near-wall  region  and  the  outer  fiow.  The  near-wall  concentration 
decreases  with  increasing  distance  downstream  until  the  concentration  is  nearly 
uniform  thirty-two  channel  heights  downstream  of  the  injectors.  This  result  is  in 
agreement  with  pressure-drop  measurements  which  indicate  that  drag-reduction  is 
nearly  constant  downstream  of  this  location. 

The  flow  has  reached  a  near  fully-developed  condition  at  the  far  downstream 
measurement  location  (x/h  =  32)  where  the  level  of  drag  reduction  was  measured  at 
thirty-four  percent.  Figures  6  and  7  compare  the  first  and  second  moments  of  the 
streamwise  and  normal  velocity  components  with  those  of  Luchik  (1985)  which  were 
made  in  the  same  channel  with  a  drag-reduction  level  of  thirty-one  percent  for  a 
fully-developed  flow  at  a  Reynolds  number  of  22,000.  The  quantities  are  normalized 
with  the  local  (drag-reducing)  values  of  shear  velocity  and  kinematic  viscosity.  With 
the  exception  of  the  root-mean-square  of  the  normal  velocity  component,  the 
agreement  between  these  two  data  sets  is  very  close.  While  these  two  flows  are  only 
similar  in  terms  of  the  drag  reduction  level,  this  comparison  indicates  that  the  flow  is 
approaching  a  condition  that  one  would  expect  in  a  fully  developed  drag-reducing 
channel  flow. 
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Figure  5 


Mean  polymer  concentration  profiles  downstream  of 
injector;  □,  x/h  =  4;  8;  A,  16;  V,  32. 


Turbulence  statistics  are  usually  presented  normalised  with  the  local  values  of 
the  shear  velocity  and  kinematic  viscosity.  For  a  Newtonian  flow  this  normalization 
results  in  "universal"  behavior  of  turbulence  statistics  in  the  near-wall  region.  Similar 
universality  has  not  been  demonstrated  for  even  fully  developed  drag-reducing  flows. 
Hence,  this  normalization  does  not  provide  a  basis  for  comparison  and  tends  to 
obscure  the  absolute  changes  in  the  statistical  quantities. 

An  absolute  measure  of  the  effect  of  the  injected  polymer  solution  on  the 
structure  of  the  flow  results  from  comparing  the  turbulence  statistics  of  the  injected 
flow  with  those  of  the  water  flow  at  the  same  physical  location.  The  following  figures 
show  the  changes  in  structure  by  presenting  the  ratio  of  a  given  quantity  for  the 
drag-reducing  flow  to  the  corresponding  quantity  for  the  same  physical  location  in  the 
water  flow  with  no  injection. 

Figure  8  shows  the  fractional  change  in  the  mean  ^treamwise  velocity,  U,  with 
polymer  injection  from  that  of  the  water  flow  as  a  function  of  distance  normal  to  the 
wall  for  the  four  streamwise  measurement  locations.  The  distance  from  the  wall,  y^, 
is  normalized  with  the  shear  velocity  and  kinematic  viscosity  of  the  water  flo'v  for 
reference  only  and  is  not  meant  to  imply  that  the  flow  structure  should  be  correlated 
by  this  normalization.  The  mean  velocity  decreased  in  the  high  polymer  concentration 
near-wall  region,  y^  less  than  one  hundred.  The  effect  is  most  pronounced  in  the 
near-injector  location  where  the  polymer  concentration  is  highest  and  diminishes  with 
increasing  streamwise  distance.  A  slight  acceleration  of  the  outer  flow  results  from 


mass  conservation  considerations. 
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The  fractional  change  in  the  root-mean-square  (RMS)  of  the  streamwise  velocity, 
u',  (Figure  9)  exhibits  similar  behavior  for  less  than  thirty  but  shows  an  increase 
for  y^  from  thirty  to  two  hundred.  The  decrease  immediately  adjacent  to  the  wall 
and  the  increase  further  from  the  wall  can  be  attributed  to  a  corresponding  decrease 
and  increase  in  the  gradient  ofU  in  these  regions.  Figure  10  shows  the  behavior  of  the 
RMS  of  the  normal  velocity,  v',  for  the  four  measurement  locations  during  polymer 
injection.  There  is  a  decrease  in  v'  across  the  entire  channel  at  all  streamwise 
locations  except  x/h=4.  At  this  location,  v'  is  reduced  only  in  the  high  concentration 
near-wall  region  and  appears  to  be  unaltered  in  the  outer  flow. 

The  fact  that  the  turbulent  shear  stress,  uv,  approaches  sero  at  the  channel 
centerline  makes  comparison  of  this  quantity  in  a  similar  manner  somewhat  more 
ambiguous.  For  this  reason,  the  comparison  is  made  by  plotting  the  measured  value  of 
uv  for  the  injected  flow  versus  distance  from  the  wall  along  with  a  line  representing  uv 
for  the  fully  developed  water  flow.  These  data  are  presented  in  Figure  11  where  both 
quantities  are  normalized  with  the  shear  velocity  of  the  water  flow.  In  a  manner 
similar  to  that  of  v',  the  value  of  uv  is  reduced  across  the  entire  channel  height  at  all 
locations  except  at  x/h  =  4.  Here  the  value  off  uv  deviates  signiBcantly  from  that  of 
the  water  flow  only  in  the  region  of  high  polymer  concentration. 

2.3  Conclusions  and  Recommendations 

From  these  data  several  general  results  can  be  identified.  The  structure  of  the 
near  wail  region,  less  than  one  hundred,  appears  to  be  highly  dependent  upon  the 
polymer  concentration.  The  mean  streamwise  velocity  gradient  is  decreased 


Figure  9.  Fractional  change  in  the  RMS  of  the  streamwise  velocity 
with  polymer  injection;  □,  x/h  =  4;  8;  A,  16;  V,  32. 
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Figure  10.  Fractional  change  in  the  RMS  of  the  normal  velocity  with 
polymer  injection;  □,  x/h  =  4;  Q,  8;  A  16;  V,  32. 
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11.  Reynolds  stress  profiles  for  flow  with  polymer  injection;  □ 
x/h  =  4;  O,  8;  A,  16;  V,  32; - ,  water  flow. 


immediately  adjacent  to  the  wall  and  increased  at  distances  further  from  the  wall. 
This  results  in  a  corresponding  decrease  and  increase  in  u*  relative  to  the  same 
location  in  the  water  flow.  The  RMS  of  the  normal  velocity  and  the  turbulent  shear 
stress  are  always  less  than  or  equal  to  the  values  for  the  water  flow.  At  x/h  =  4, 
these  quantities  deviate  greatly  from  the  corresponding  water  flow  levels  only  in  the 
high-concentration  layer.  Further  downstream,  where  the  polymer  is  present  across 
the  entire  channel  height,  V  and  uv  are  reduced  everywhere.  This  indicates  that  once 
the  polymer  has  mixed  beyond  the  wall  region  in  even  small  quantities  the  structure  of 
the  flow  has  been  altered  in  every  aspect. 

Since  the  spanwise  average  of  the  dye  concentration  across  the  channel  was  not 
constant  with  time,  the  attenuation  method  cannot  be  used  to  deduce  instantaneous 
values  of  the  polymer  concentration  (see  Tiederman  et  al.,  1986).  Instead,  time- 
resolved  concentration  measurements  will  be  made  using  the  fluorescent  dye  technique 
of  Koochesfahani  and  Dimotakis  (1986).  One  objective  of  our  future  efforts  will  be  to 
measure  directly  the  product  of  the  normal  velocity  fluctuation,  v  and  the  polymer 
concentration  fluctuation,  c^  The  time  average  of  this  product  is  the  additional 
unknown  in  the  Reynold’s  averaged  equation  for  polymer  concentration. 

The  sharp  gradients  of  the  average  polymer  concentration  imply  that  even 
sharper  gradients  exist  instantaneously.  Future  efforts  should  focus  on  this  high 
Schmidt  number  mixing  in  the  near-injector  region. 


3.  NEWTONIAN  BURST  FREQUENCY 


The  objective  of  this  experiment  was  to  confirm  the  extrapolation  made  by 
Luchik  and  Tiederman  (1987)  that  the  average  time  between  bursts,  Tq  (the  inverse  of 
the  burst  frequency)  scales  with  inner  variables  for  Newtonian,  fully  developed  channel 
fiows.  Specifically  the  proposed  relationship  is 


Tb 

v 


90 


for  Reynolds  numbers  >  10,000. 


(2) 


3.1  Apparatus  and  Procedures 

The  experiments  were  conducted  in  the  2.5  X  25.0  cm  channel  and  fiow  loop 
described  in  Section  2.1.  No  fiuid  was  injected  through  the  slots  and  the  water  was 
circulated  at  the  maximum  fiow-rate  capacity  using  all  four  pumps.  For  this  situation 
the  Reynolds  number  based  on  the  average  fiow  rate  and  channel  height,  h=2.5  cm, 
was  nominally  50,000. 

Streamwise  velocity  measurements  were  made  at  y'*'=30  by  operating  the  laser 
velocimeter  on  the  green  line  (514.5  nm)  in  a  one-component,  dual  beam  configuration. 
This  one-component  configuration  was  chosen  because  it  was  clear  prior  to  the 
experiments  that  it  would  be  difficult  to  achieve  data  rates  equal  to  or  greater  than 
the  viscous  frequency  {xx^ju  =  8,800  Hz)  with  a  spanwise  probe  volume  length  z'''<20. 
These  constraints  on  data  rate  and  probe  volume  length  were  suggested  previously  by 
Bogard  and  Tiederman  (1986)  and  Luchik  and  Tiederman  (1985,  1987).  Since  all  high 


Reynolds  number  laser  velocimeter  measurements  are  limited  by  these  constraints, 
data  were  also  acquired  with  a  larger  probe  volume  and  a  slower  sampling  frequency 
in  an  effort  to  define  these  constraints  more  precisely. 

All  of  the  water  in  the  flow  loop  was  seeded  uniformly  with  2.02  /im  diameter, 
white,  latex  spheres.  These  scattering  particles  were  chosen  because  their  specific 
gravity  is  nearly  1.0  and  the  fringe  spacing  of  the  laser  velocimeter  was  2.68  /im. 

The  receiving  optics  were  aligned  to  receive  forward  scattered  light  at  an  angle 
of  either  0  *  or  45  ‘  with  respect  to  the  optical  axis  of  the  transmitting  optics.  At  45  * , 
the  probe  volume  length  was  truncated  to  a  spanwise  length  given  by  z'*'  =21.  At  0  * , 
the  probe  volume  length  was  z"*"  =  137. 

Measurements  were  made  with  the  45*  configuration  at  particle  number 
densities  that  produced  nominal  data  rates  of  4  kHz  and  8  kHz.  The  receiving  optics 
were  then  realigned  to  0*  and  the  particle  density  was  adjusted  to  produce  an  8  kHz 
data  rate.  After  collecting  and  storing  data  from  the  free-running  counter  processor, 
the  velocity  records  were  then  sampled  at  equal  intervals  of  time  corresponding  to 
their  respective  data  collection  rates  using  linear  interpolation  between  data  points. 

The  techniques  for  grouping  single-component  velocity  detections  of  ejections 
into  burst  detections  as  presented  by  Luchik  and  Tiederman  (1987)  were  used  in  this 
study.  In  particular  the  u-level,  modified  u-level  and  VITA  algorithms  were  used  to 
detect  ejections.  The  u-level  method  has  one  threshold  parameter,  L,  and  a  detection 
occurs  when  the  fluctuation,  u,  about  the  mean  U  is  given  by 
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u  <  — L  u'  (3) 

Here  u'  is  the  root-mean-square  of  u.  For  the  modified  u-level  technique  Equation  3 
describes  the  condition  when  the  detector  is  turned  "on."  A  modified  u-level  detection 
ends  when 


u  >  —0.25  L  u' 

The  threshold  for  a  VITA  detection  is  given  by 

VAR  >  k  u'* 

where 

VAR  =  0^-0* 

and 


(4) 


(5) 


(6) 


t+'-<.TA 

/  Udt  (7) 

The  value  for  Ty^  was  chosen  to  yield  the  maximum  number  of  detections  independent 
of  threshold.  Ejections  are  then  grouped  into  burst  detections  by  using  the  probability 
functions  for  time  between  ejections  to  estimate  a  time  period,  that  separates 
ejections  from  the  same  burst  from  ejections  from  other  bursts.  If  the  time  between 
temporally  adjacent  ejections  is  less  than  rg,  those  two  ejections  are  considered  to  be 
part  of  one  burst.  Conversely  if  the  time  between  adjacent  ejections  is  larger  than  fg, 
then  those  two  ejections  are  from  different  bursts. 


3.2  Results 


I  - 

1 

Figure  12  which  is  a  plot  of  the  probability  that  T  >  Tg  demonstrates  the 
procedure  for  determining  for  u-level  detections.  This  is  the  same  technique  used  by 
j  Barlow  and  Johnston  (1985)  and  Luchik  and  Tiederman  (1987)  for  quadrant  2 

detections  from  two-component  data.  As  indicated  by  the  different  symbols  on  Figure 
13,  the  value  of  rg  varied  somewhat  with  specification  of  the  threshold  L.  The  u-level 
method  yields  an  average  time  between  bursts,  Tg,  of  8.0  ms  for  the  data  set  where 
z'*'  =  21  and  the  sampling  data  rate  equals  8  kHz. 

The  u-level  and  quadrant  2  procedure  for  determining  rg  does  not  work  for  either 
modified  u-level  or  VITA  detections.  The  procedure  of  leaving  the  detector  "on"  in  the 
modified  u-level  method  and  the  integration  time,  T^,  in  the  VITA  method  effectively 
low-pass  filter  the  velocity  record.  As  a  result  both  methods,  by  design,  do  not  detect 
all  of  the  ejections  from  the  same  burst.  Consequently  the  probability  function  is 
reduced  substantially  for  small  values  of  time  between  ejections,  Tg,  while  it  is 
unaffected  at  large  values  of  Tg.  Consequently,  a  method  based  on  fitting  the 
histogram  of  the  number  of  ejections,  Ng,  in  a  small  interval  of  Tg  to  a  Poisson 
distribution  was  used  to  determine  rg. 
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Figure  12.  Probability  function  for  time  between  ejections  using 
u-level  detection;  t*  =  21,  8  kHz  data. 
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Figure  13.  Time  between  bursts  for  u-level  method; 

O.  'e  — 2  5  '■£  =  2.7  ms;  V.  "e  =  2.6  ms; 

2.4  ms;  ,  '^e^2.1  ms. 
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In  practice  this  procedure  for  matching  the  experimental  histogram,  shown  in 
Figure  14,  to  the  theoretical  distribution  is  dependent  upon  the  investigators  choice  of 
"bin"  site  which  in  this  case  is  the  iutervsl  in  Tg.  If  the  interval  is  too  small  the 
histogram  does  not  approximate  a  smooth  curve  while  if  the  interval  is  too  large  the 
method  has  poor  resolution.  With  several  hundred  ejection  detections,  a  good 
compromise  is  to  choose  a  time  interval  which  places  the  mode  in  the  fourth  or  fifth 
bin.  The  mode  of  this  experimental  distribution  is  then  set  equal  to  the  mean  value  of 
the  Poisson  distribution  and  Cg  is  chosen  to  be  the  value  of  Tg  where  P(T<Tg)  *  0.85 
for  the  Poisson  distribution. 

An  illustration  of  the  grouping  technique  for  modified  u-level  deteciuins  is  shown 
in  Figure  15.  This  short  segment  of  the  data  record  was  chosen  because  it  contains 
examples  where  the  method  worked  well  as  well  as  examples  where  it  may  not  have 
properly  grouped  detections.  For  example,  based  on  the  u  signal,  the  group  of  3 
ejections  between  170  and  175  ms  are  most  likely  one  burst.  On  the  other  hand  the 
two  ejections  between  156  and  160  ms  are  most  likely  from  different  bursts  because  the 
streamwise  velocity  is  so  much  larger  than  the  mean  value  in  the  interval  between  the 
detections.  We  did  not  expect  these  simple  grouping  methods  to  be  perfect. 
Undoubtedly  additional  criteria  could  be  developed  that  would  improve  the  precision 
of  the  grouping.  At  present  the  uncertainty  that  propagates  to  the  determination  of 
Tg  is  probably  on  the  order  of  20%.  Despite  this  uncertainty,  which  is  common  to  all 
statistical  methods  for  estimating  Cg,  the  modified  u-level  method  yields  a  value  of 
Tg  =  12.0  ms  over  a  rather  wide  range  of  threshold  levels  for  the  i'  =  21,  8  kHi  data 
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The  Poisson  method  for  estimating  'g  must  also  be  used  when  the  ejections  are 
detected  using  the  VITA  method.  As  shown  in  Figure  17,  the  procedure  yields 
Tg  *  10.0  over  a  small  but  finite  range  of  the  threshold  parameter,  k,  for  the  i*  =21, 
8  kHi  data. 

The  average  value  of  Ty  from  the  three  different  methods  used  to  analyze  the 
1*  =21.  8  kHz  data  was  10.0  rns.  The  same  three  methods  were  applied  to  the 
X*  *21,  4  kHz  data  set  and  the  i*  =  137,  8  kHz  data  set.  It  was  more  difficult  to 
determine  for  these  data  sets  l>ccause  some  temporal  resolution  was  missing  from 
the  4  kHz  data  and  some  spatial  averaging  occurred  in  the  z^  ^  137  data.  These 
"degradations  "  mean  that  some  ejections  are  not  detected  by  the  procedures  and  this 
affects  plots  such  as  those  on  Figures  12  and  14.  Nonetheless  estimates  were  made 
using  all  three  methods  and  the  average  value  of  Tg  for  these  two  additional  data  sets 
are  reported  in  the  nondimensions  I  plots  that  show  the  variation  of  Tg  with  Reynolds 
number. 

Figures  18,  19,  and  20  show  how  the  present  data  along  with  previous  data 
(Luchik  and  Tiederman,  1987)  scale  with  Reynolds  number.  All  of  these  results  were 
obtained  in  the  2.5  x  25.0  cm  channel.  The  outer  scaling  of  Tg  using  the  renter-line 
velocity,  Uq,  and  channel  half-height,  h/2.  clearly  does  not  yield  a  constant 
nondimensional  value.  Neither  does  the  mixed  scaling  proposed  by  Alfredsson  and 
Johanssen  (1984).  Scaling  Tg  with  inner  variables  does  work  yielding 

fg  =  90 


3.3  Conclusions 


The  principal  conclusion  is  that  T3  =  90  for  9,400  <  Rej,  <  50,000  in  a  fully 
developed  channel  flow.  Based  on  our  experience  with  data  from  longer  spanwise 
probe  volumes  and  with  data  rates  less  than  the  viscous  frequency,  it  is  recommended 
that  the  probe  volume  have  a  spanwise  length  of  z'''  <  20  and  that  the  data  rate  be 
equal  or  greater  than  the  viscous  frequency,  Ur/i^. 
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Figure  18.  Outer  scaling  for  time  between  bursts;  Oi  Luchik  & 
Tiederman  (1987);  ,  z+=21,  4  kHz;  A,  z‘"=21,  8 

kHz;  □,  z'*'=137,  8  kHz. 
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Figure  19.  Mixed  scaling  of  time  between  bursts  (see  Figure  18 
for  symbols). 
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